Gradient Echo Quantum Memory for Light using Two-level Atoms 
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We propose a quantum memory for light that is analogous to the NMR gradient echo. Our 
proposal is ideally perfectly efficient and provides simplifications to current 3-level quantum memory 
schemes based on controlled inhomogeneous broadening. Our scheme does not require auxiliary 
light fields. Instead the input optical pulse interacts only with two-level atoms that have linearly 
increasing Stark shifts. By simply reversing the sign of the atomic Stark shifts, the pulse is retrieved 
in the forward direction. We present analytical, numerical and experimental results of this scheme. 
We report experimental efficiencies of up to 15% and suggest simple realizable improvements to 
significantly increase the efficiency. 
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PACS numbers: 

Some of the most significant advances in quantum in- 
formation processing have been made using quantum op- 
tics. To extend these techniques, it is necessary to have 
devices such as single photon sources, quantum memo- 
ries and quantum repeaters, where quantum information 
is exchanged in a controlled fashion between light fields 
and material systems. 

A quantum memory for light is a device that can effi- 
ciently delay or store the quantum states of light fields. 
This is usually achieved via some form of imprinting onto 
an atomic system. The quantum states stored by the de- 
vice must also be faithfully retrievable on demand where 
the total efficiency of the processes must exceed the clas- 
sical benchmark so that quantum information can be re- 
tained [lj. It has been proposed that the requisite con- 
trol and strong coupling can both be achieved using an 
ensemble approach, where the light field interacts with 
a large number of identical atoms [2]. In the discrete 
variable regime, single photon states have been stored 
and retrieved from magneto-optical traps using electro- 
magnetically induced transparency [U [5] . Whilst in the 
continuous variable regime, quantum states of light have 
been mapped onto atoms using the resonant interaction 
of light with spin polarised cesium vapors [3]. Never- 
theless, further improvements on the efficiencies are still 
needed for realizing a reliable quantum memory. 

In 2001, Moiseev and KrdLl [6 published a proposal 
for a quantum memory for light based on modified pho- 
ton echoes. In contrast to a normal photon echo, the 
rcphasing came from controlled reversible inhomogeneous 
broadening (CRIB) rather than from optical pulses. In 
the initial proposal the reversible inhomogeneous broad- 
ening used opposite propagation directions in a Dopplcr 
broadened medium. Since then the idea has been gener- 
alised to other broadening mechanisms [Tj |5J [5] ■ 

Photon echoes using reversible inhomogeneous broad- 
ening [10 , with the ability to store multiple pulses, 
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FIG. 1: (Color online) a) Schematic of a gradient echo quan- 
tum memory, (i) The input light field enters the sample at 
z — —Zo and t — —to; (ii) the optical information of the 
light field is imprinted onto the inhomogenously broadened 
two level atomic ensemble, (iii) to recall the light field, the 
polarity of the quadrupole electric field is flipped at t = 0, 
(iv) finally an output field is retrieved at z — zo and t = to- 
With the optical depth chosen to be gN = 3.3?? and with 
no decoherence 7 = 0, b) and c) show the space-time grid 
plots of the light field intensity and the atomic polarization, 
respectively. The input pulse duration is t pu i sc = to /A and 
the quadrupole induced broadening is 2/t pu i sc . 



has been demonstrated using Stark shifts in europium 
dopants [TTJ. Work towards demonstrating such echoes 
in other systems has also been carried out [SJ [T2J [T3] . To 
date all the proposals for a quantum memory using CRIB 
operate via a time reversal of the storage process. By re- 
versing the detunings of the atoms, one can transform 
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the equations of motion for light travelling in the back- 
ward direction into a time reversed copy travelling in the 
forward direction. A forward propagating input pulse is 
absorbed by the atomic ensemble. Then the detunings of 
the atoms are flipped and a phase matching operation is 
applied. The phase matching operation consists of a pair 
of counter propagating it pulses driving the atoms from 
the excited state down to, and back up from, an auxiliary 
third level. Finally the pulse exits the ensemble in the 
backward direction as a time reversed copy of the input. 

In this paper, we show that 100% efficiency using CRIB 
is theoretically possible using only two level atoms. Fur- 
thermore the output pulse propagates in the same direc- 
tion as the input pulse. It had been thought that such 
two-level schemes would not be able to approach unit ef- 
ficiency because of the problem with re-absorption. We 
will show, however, that this is not the case provided 
that the inhomogeneous broadening is introduced in a 
way that the detuning of the atoms is linear with posi- 
tion. The principle benefit of such a two- level scheme is 
in its simplicity. Firstly, the absence of phase matching 
7r pulses greatly simplifies the implementation. The only 
light seen by the atomic ensemble during the operation 
of the memory is the light field of interest. Secondly, 
the memory requires only two atomic levels making this 
scheme applicable to many more atomic systems. In par- 
ticular erbium dopants which allow operation at the tele- 
communication wavelength of 1.5 /im, have been shown 
to have very good two-level characteristics [14], whilst a 
lambda system has yet to be demonstrated. 




FIG. 2: (Color online) Contour plot of the real part of the 
optical field in a moving frame at the speed of light. From 
t — —to to t = 0, the light field enters the sample and is 
absorbed by the medium. At t = 0, the quadrupole field is 
flipped and the time reverse process commences producing a 
forward propagating pulse. For the same parameters as given 
in Fig. [I] with a storage time 8 x the pulse duration (i) shows 
a small phase shift across the retrieved pulse. With a storage 
time 80 x the pulse duration, (ii) shows a near ideal pulse 
retrieval. The arrows denote the wavefront of the light fields. 

We consider the interaction of a collection of two- 
level atoms with a light field where the detuning of the 



atoms is linearly dependent on their position. After the 
small pulse approximation, which ensures that a neg- 
ligible amount of atoms reaches the excited state, the 
Maxwell-Bloch equations in the frame at the speed of 
light are [TH] 
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Where E represents the slowly varying envelope of the 
optical field; a the polarisation of the atoms; N the 
atomic density; g the atomic transition coupling strength; 
7 the decay rate from the excited state and rjz the de- 
tuning from resonance. Although the treatment here is 
classical, the linearity of Eqs. (Tj2l ensures that the re- 
sults will hold for quantised fields. 

As shown in Fig.[l]a), the pulse enters the medium at 
z = —zq, t = —to- The detunings of the atoms is flipped 
(—irjz — > +irjz) at t = 0. We will show that this leads to 
a forward propagating output pulse that is closely related 
to the input at z — z , t — t Q . Fig. [l]b) shows the results 
of numerical simulations for the gradient echo scheme. 
For these simulations and the analytical discussion below, 
we assumed a zero decay rate 7 and a large optical depth 
gN/rj = 10/3 with an inhomogeneous broadening two 
times larger than the spectral width of the input pulse. 
In practice such a narrow linewidth two-level atom can 
be realized by optical pumping to a third level that does 
partake in the quantum memory process. Under these 
conditions, the spectral components of the pulse are ab- 
sorbed at different positions in the sample, which results 
in a temporal stretching of the pulse as it propagates. 
When the sign of the inhomogeneous broadening is re- 
versed, the atoms re-radiate the individual spectral com- 
ponents as they re-phase. Because of the monotonicity of 
the Stark shift, the time reversed copy of the input pulse 
propagate in the forward direction with no re-absorption. 

Fig. [2] shows a contour map of the real part of the 
electric field. Because of the large phase shift seen by 
the field, the last atoms in the sample absorb the field a 
long time after the first atoms. When the switching time 
is short relative to the pulse duration, a time dependent 
phase shift is present across the output pulse causing it 
to be frequency shifted with respect to the input pulse. 
Fig. [2ji) shows results where the storage time is eight 
times the pulse duration. For long switching time, nearly 
all the light field is imparted onto the atomic ensemble. 
This leads to a uniform phase shift for the entire pulse 
and the production of a replica of the input field as shown 
in Fig. gii). 

To gain a better understanding of the process we 
present an analytic treatment which gives an expression 
for the output field as a function of the input i n the case 
of large optical depth. First we solve solve Eqs (l]2l sub- 
ject to the boundary conditions of a(z, t — > —00) = 
and E(z — —z ,t < 0) — f m (t) to find the values for 
a and E in the region (t S (— oo,0], z e [— zq,zq]). To 
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find the values for a and E in the region (t S [0, oo), 
z G l—Zp , zp]) one needs to solve a modified version of 
Eqs (Tj2| where the sign of the ir/z term is reversed. One 
way to do this would be to propagate these equations for- 
ward. However we will take a different approach which 
makes better use of the symmetry of the protocol. We 
will guess a particular form of the final conditions, that 
is E(z = ZQ,t > 0) = /out(t), a(z,t — > oo) = and prop- 
agate these expressions backwards to find a(z,t = 0), 
E(z,t — 0). We will then validate that our guess is cor- 
rect by showing that the boundary conditions for our 
backward propagated solutions for the region t > match 
those that we found for the region t < 0. Integrating 
Eq. ([!]) gives 



a(z,t) =ig dt' H(t-t')e 



— ir)z(t—t ) 



E{z,t') (3) 



where H(£) is the Heaviside step function. Fourier trans- 
forming this expression, substituting into Eq.Q and then 
integrating the result gives 



E(z,u) = F in (u) exp[-TT/3(H(u + -qz) - H{w - t]Z )) 
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FIG. 3: Energy level diagram and spectral scheme for the 
praseodymium dopants in yttrium orthosilicate. The gradient 
echo is carried out with light at frequency A. The natural 
inhomogeneous linewidth of the sample, w n at, is a few GHz 
wide. To set up the experiment, the applied light is swept 
around frequency A to create a spectral hole a few MHz wide, 
Whole- Light at frequency B and C is then applied to prepare 
a narrow antihole around A with linewidth w an tj. Although 
the diagram shows the use of the ±1/2 excited states, there 
are ions contributing to the antihole from the ±3/2 and ±5/2 
excited states due to inhomogeneous broadening in the optical 
transition. In our experiment u> an tihoie = 30 kHz. 



Here Fi n (uS) is the Fourier transform of fi n (t) and (3 — 
gN/r] is the optical depth. It can be seen that each spec- 
tral component is attenuated by a factor exp(— (3) after 
travelling past the position in the sample where it is res- 
onant with the atoms, as well as getting a phase shift as 
it travels through the sample. 

Because of this phase shift, the intensity of the field 
tends to zero with time asymptotically. Its expression at 
t = 0, can be obtained by integrating Eq.Q over all u>, 
where we assume that the bandwidth of the memory is 
larger than the bandwidth of the input pulse, t/zq u>. 
This simplification leads to an expression for E{z,t = 0) 
in the form of a convolution. Fourier transforming with 
respect to the spatial coordinate leads to 
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Here fin{k/rf) is the input field at the time r = k/rj and 
r(£) is the gamma function. 

As mentioned above we can also derive an expres- 
sion for E(k, t = 0) subject to the final conditions of 
a(z,t — > oo) = and E{z — zo,t > 0) = f ou t(t), 
one gets a similar expression to Eq.(|5| but in terms of 
fout(k/i]). The backwardly propagated solution satis- 
fies the boundary condition E(z = —Zo,t > 0) = and 
a(z,t — > oo) = so long as the optical depth is suffi- 
ciently high. Using continuity arguments, we match the 
two solutions found for E{k, t = 0) and get an expression 
of the output field at a time to. 



/out ^o) — /in(-to)|*o| 2l/3 -p/ T5\ ■ 

T{-i(3) 
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Because \to\ 2tf3 T(i(3)/T(— i(3) has a modulus of 1 for 
all t , it can bee seen that the envelope of the output 
pulse is a time reversed version of the input. When the 
length of the pulse is small compared to the total delay 
it is also phase shifted by |to| • This term describes a 
frequency shift of the output pulse with respect to the 
input. One other elegant way to cancel this frequency 
shift is to cascade two gradient echo memories where the 
flipping is performed the opposite way for the second 
memory. It can be see from Eq.(6]) that the frequency 
shift from each memory would indeed cancel. 

The theoretical treatment presented above showed that 
this Gradient Echo scheme can in principle be used as 
a quantum memory for light. Preliminary experiments 
have also been carried out in a similar manner to the 
initial demonstrations of photon echoes produced using 
controlled inhomogeneous broadening |10j . The most sig- 
nificant difference was the use of praseodymium rather 
than europium dopants, which provided much greater 
optical depths. The experiment was carried out on a 
spectral antihole which was prepared as described in 
Fig. [3j Light from a highly stabilised dye laser was fre- 
quency shifted and gated with acousto-optic modulators 
(AOMs). The pulse was then steered toward the sample 
of 0.05 atomic % praseodymium doped in yttrium or- 
thosilicate. The sample was approximately a 4 mm cube 
and was held at temperatures in the range 2-4 K. Four 
electrodes were placed around the sample in a quadrupolc 
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FIG. 4: Experimental results (solid lines) and the numerical 
simulations (dash lines) of the gradient echo setup. When 
no antihole is prepared, trace (a) shows the input pulse being 
detected after transmission through the sample. The gradient 
echo trace shows (b) the transmitted and (c) the stored-and- 
recalled pulse components. The vertical line at 3.7 fis denotes 
the time at which the Stark shifting electric field is flipped. 

arrangement and provided an electric field that varied lin- 
early along the optical path. The electrodes were 1.7 mm 
diameter rods separated by 8 mm. Voltages of approxi- 
mately ±5 V were used to broaden the antihole. Hetero- 
dyne detection was used to detect the transmitted pulses. 

Fig. [4] shows the typical experimental traces of the gra- 
dient echo memory with and without the preparation of 
the two-level anti-hole. In the results, 49% of the inci- 
dent light was transmitted straight through the sample 
and 13% of light was recalled as an echo. Due to some 
variation in the antiholes created, we observed a variation 
of the efficiency of the echo between the values of 10 and 
15%. Using our numerical model, we vary the spectral 
width of the unbroadened anti-hole and the optical depth 
of the sample to match the experimental results. Close 
agreement between the experimental results and the sim- 
ulations is obtained only with these two free parame- 
ters. Our numerical model suggests an anti-hole width of 
30 kHz. This is again in agreement with the experimen- 



tal expectation, if the hyperfine transition broadening 
were the main limitation to the antihole linewidth. For 
a given pulse length, the optimization of the experiment 
is dependent on a compromise between maximizing the 
optical depth of the sample and increasing the linewidth 
ratio between the applied and the intrinsic inhomogc- 
nous broadenings. In our experiment, gN/r/ = 0.006, 
whilst the ratio between applied and intrinsic linewidth 
is around 200. 

For each crystallographic site where praseodymium is 
located, there is another related to it by inversion. In or- 
der to implement a completely efficient memory only one 
of the site pair can be used. In principle, this could be 
achieved by Stark shifting with a homogeneous electric 
field and optical pumping. In our experiment, however, 
both orientations were used. The theoretical modeling on 
Fig.[4]takes into account these two orientations by having 
two Bloch equations and two source terms for the opti- 
cal field. The performance of the gradient echo memory 
can be enhanced with minor improvement to the set-up. 
Simulations suggest that with Fourier limited pulses, se- 
lecting only one orientation of the praseodymium ions 
and increasing the optical depth by a factor of three (for 
instance, by increasing the crystal length) would enable 
the scheme to reach more than 50% efficiency. 

In conclusion we have proposed a quantum memory 
for light based on optical gradient echoes. When com- 
pared to existing quantum memories based on controlled 
inhomogeneous broadening, our scheme requires only two 
atomic levels and is therefore easier to implement on 
a wider range of systems. Moreover, our scheme does 
not require optical pulses for the storage and recall pro- 
cess. Initial experiments show an average echo efficiency 
of 13% and a time bandwidth product of around three. 
This compares favourably with the performance of quan- 
tum memories based on EIT [H [5] . The experiments are 
well modeled by simple Maxwell-Bloch equations. Simu- 
lations of modest improvements on the experimental pa- 
rameters suggest that efficiency higher than 50% should 
be realizable. 
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